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ABSTRACT 

Differentially rotating, 'advection-dominated' accretion flows are considered in which 
the heat generated by viscous dissipation is retained in the fluid. The equations of 
time-dependent quasi-spherical accretion are solved in a simplified one-dimensional 
model that neglects the latitudinal dependence of the flow. A self-similar solution is 
presented that has finite size, mass, angular momentum and energy. This may be 
expected to be an attractor for the initial-value problem in which a cool and narrow 
ring of fluid orbiting around a central mass heats up, spreads radially and is accreted. 
The solution provides some insight into the dynamics of quasi-spherical accretion and 
avoids many of the strictures of the steady self-similar solution of Narayan & Yi. 
Special attention is given to the astrophysically important case in which the adiabatic 
exponent 7 = 5/3; even in this case, the flow is found to be differentially rotating and 
bound to the central object, and accretion can occur without the need for powerful 
outflows. 
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1 INTRODUCTION 

Accretion flows in which a significant fraction of the heat 
generated by viscous dissipation is retained in the fluid, 
rather than radiated away, have been the subject of con- 
siderable attention in recent years (e.g. Narayan & Yi 
1994, 1995a,b; Abramowicz et al. 1995). These 'advection- 
dominated' flows occupy an intermediate position between 
" the spherically symmetric accretion flow of a non-rotating 
fluid (Bondi 1952) and the cool, thin disc of classical 
accretion-disc theory (e.g. Pringle 1981). Although their self- 
consistency may be contingent on certain unresolved issues 
of plasma physics, they have been considered relevant to 
a variety of astrophysical objects, in particular to explain 
the low luminosity of the Galactic Centre (e.g. Narayan et 
al. 1998). A simple algebraic model describing the fluid- 
dynamical aspects of such a flow, which was assumed to be 
steady, axisymmetric and radially self-similar, was presented 
by Narayan & Yi (1994, hereafter NY). They adopted a set 
of approximate spherically (or vertically) averaged equations 
which were largely justified by subsequent calculations that 
included the latitudinal dependence (Narayan & Yi 1995a; 
see also Gilham 1981). 

Two properties of this solution have been emphasized 
and identified as potential difficulties with its application. 
First, in the astrophysically important case in which the 
adiabatic exponent 7 = 5/3, the solution of NY coincides 
with that of Bondi (1952). This is spherically symmetric 
and non-rotating, has no viscous torque or dissipation, and 
has been considered to depart too far from the concept of an 



accretion disc. In particular, it is not clear how such a solu- 
tion could connect to the initial or external conditions if the 
supplied matter had any angular momentum. To circumvent 
this difficulty, some authors have investigated the possibil- 
ity that the effective value of 7 might be less than 5/3 if the 
magnetic field contributed a significant fraction of the total 
pressure (Narayan & Yi 1995b; Quataert & Narayan 1998). 

Secondly, the Bernoulli parameter of the flow is posi- 
tive. This has been interpreted as meaning that the fluid is 
unbound and would generate a powerful outflow if the con- 
straints of the model were somehow relaxed. Xu & Chen 
(1997) therefore generalized the work of Narayan & Yi 
(1995a) by constructing self-similar solutions that involve 
inflow at equatorial latitudes and outflow at polar latitudes. 
However, to achieve this, the scalings of physical quanti- 
ties with radius were adjusted arbitrarily from their natural 
values so that only an asymptotically small fraction of the 
mass reaches the central object, and in this sense the flow is 
hardly an accretion flow. The inflow is forced to turn into an 
outflow because the central object by assumption accretes 
essentially no mass. These solutions are relatively uncon- 
strained because they are unbounded in size and energy and 
are not required to satisfy any radial boundary conditions. 
It is therefore uncertain whether they could be realized as 
intermediate asymptotic forms in practice. 

In a recent paper, Blandford & Begelman (1999) have 
sought to resolve both these difficulties by elaborating on the 
basic model, adding terms that parametrize the removal of 
mass, angular momentum and energy by an outflow, which 
is considered necessary 'to allow accretion to proceed'. In 
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their solution the flow is rotating even though 7 = 5/3, and 
the Bernoulli parameter can be negative. However, this was 
again achieved by adjusting the scalings so that essentially 
no mass reaches the central object. 

The purpose of this letter is to investigate the fluid 
dynamics of quasi-spherical accretion from an alternative 
viewpoint which avoids many of the strictures of steady self- 
similar solutions. In the theory of thin discs, a consideration 
of the initial- value problem, in which a cool and narrow ring 
of viscous fluid orbiting around a central mass spreads radi- 
ally and is accreted, provides an excellent demonstration of 
the dynamics of disc accretion and is in some ways more in- 
formative than steady solutions. In that problem the surface 
density satisfies a diffusion equation which, depending on 
the properties of the viscosity, is either linear, in which case 
the Green function can be obtained analytically (Liist 1952; 
Lynden-Bell & Pringle 1974), or non-linear, in which case 
similarity methods can be used to obtain an exact special 
solution of the equation (Pringle 1974; Lin & Bodenheimer 
1982; Lin & Pringle 1987; Lyubarskii & Shakura 1987). This 
self-similar solution accurately describes the behaviour of all 
solutions of the initial-value problem long after the inner ra- 
dius of the fluid first reaches the central object. In the case of 
quasi-spherical accretion, the governing equations are non- 
linear and similarity methods can again be used to provide 
an exact solution which may similarly be expected to be an 
attractor for the initial-value problem. 



2 BASIC EQUATIONS 

A simplified one-dimensional model for non-relativistic 
quasi-spherical accretion is adopted in which physical quan- 
tities depend only on the spherical radius r and on time t. 
The model is described by the equations 
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Here p is the density, u the radial velocity, Q the angular 
velocity, $ = —GM/r the gravitational potential (where M 
is the mass of the central object), p the pressure, £ = r 2 Vl 
the specific angular momentum, /1 the (dynamic) viscosity 
(which, for simplicity, is assumed to act only on the dif- 
ferential rotation) and e the specific internal energy. The 
Lagrangian time derivative is 
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Dt ~ dt + U dr' 
It is further assumed that the fluid is an ideal gas with 

p = (7 - l)pe, 1< 7 < 5/3, (6) 

and that the viscosity is given, as in NY, by an a-model, 



where Qk = (GM/r 3 ) 1 / 2 is the Keplerian angular velocity. 
The parameters 7 and a are assumed to be constant. 

These equations may be obtained from the full axisym- 
metric equations in spherical polar coordinates (r, 9, by 
considering the equatorial plane 9 — n/2 and neglecting 
terms associated with the 9- velocity and any ^-dependence. 
The variables may be considered as approximate spherically 
averaged quantities. Although approximate, the equations 
have two important properties: they reduce exactly to those 
of spherical accretion in the absence of rotation, and they 
possess exact conservation laws for mass, angular momen- 
tum and energy, which are of the form 
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the total energy flux. These equations are appropriate to 
quasi-spherical accretion, rather than disc accretion, because 
the dissipated energy is retained in the fluid and may be 
advected with the flow. 



3 SELF-SIMILAR SOLUTION 
3.1 Analysis 

NY solved essentially equations (0)-(^) for the case of a 
steady, radially self-similar flow. Such a solution is un- 
bounded in size, mass, angular momentum and energy, and 
must therefore be interpreted as an approximation to a fi- 
nite steady solution within a region far removed from its 
boundaries. In contrast, let us seek a time-dependent solu- 
tion of finite size, for which the integrals representing the 
total mass, angular momentum and energy, 
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are all convergent. The central object is considered to be 
arbitrarily small and to act as a sink for mass and (negative) 
energy but not as a source for angular momentum. Therefore 
J is strictly conserved, but M and E are not. 

The time-dependent self-similar solution is found by a 
well-known procedure (e.g. Barenblatt 1979). In the stan- 
dard terminology the solution is of type I, which means that 
the similarity variable can be deduced simply from dimen- 
sional considerations. In this case the dimensional constants 
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present are the quantity GM and the total angular momen- 
tum J. One therefore identifies the similarity variable 



where 



£ = r(GMt 2 )~ 1/3 

and seeks a solution of the form 

p(r,t) = P ^)J(GM)- s/s t- 7/3 , 
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In such a solution each physical quantity retains a simi- 
lar spatial form as the flow evolves, but the characteristic 
length-scale of the flow increases proportionally to t 2 ^ 3 . Note 
that the total mass M is proportional to t -1 ^ 3 and the total 
energy E (which is in fact negative) proportional to t . 

Substitution into equations (|l])-(§) yields the dimen- 
sionless equations 
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where a prime denotes differentiation with respect to £. The 
constraint 
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provides a normalization condition for the density. 

This is a fifth-order system of non-linear ordinary dif- 
ferential equations. Critical points occur wherever 
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i.e. wherever the radial velocity measured with respect to a 
self-similarly expanding coordinate system is equal to zero 
or to the local sound speed. 

3.2 Inner limit 

When 7 < 5/3, an appropriate asymptotic solution as £ — > 
is of the form 
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while A p is determined subsequently from the density nor- 
malization. The parameter A (such that < A < 1) and the 
vector [Bp, B u , Bsi, B e ] T satisfy a certain algebraic eigen- 
value problem. At leading order, we have essentially the so- 
lution of NY. This is reasonable because, for small £, the flow 
has experienced many orbits and, except for the declining 
density, may be expected to approach a steady state. 
When 7 = 5/3, we have instead 
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The limit 7 — » 5/3 is a singular one because A — > 
so that An 7^ limj4 n and A e 7^ limA e . In other words, 
the NY solution is subject to fractional corrections which 
become of order unity at all radii as 7 — > 5/3, and the Bondi 
solution is never attained. The further development of the 
inner solution consists in general of an irregular power series 
which is beyond the scope of this paper. 

3.3 Outer limit 

Any solution of finite size must possess a free surface at 
some point £ = £ s . At such a point the velocity must equal 
u — (dr/dt)j, which implies it* = §£ s . Assuming that the 
entropy tends to a finite limit, the surface resembles that of 
a polytrope. In terms of s — £ s — £, we have 

~ s 1 ^- 1 \C p +D p s + -^ (41) 

M0 ~ §& +fu«+-, (42) 

~ Cn + J Dns + --- I (43) 

e«(0 ~ s(C e + D e s + ---). (44) 

All the coefficients in this expansion, except for the density 
normalization, follow algebraically from a knowledge of £ s . 
The expansion fails if 7 < 9/7. 

3.4 Numerical solution 

If the value of £ s is guessed, the equations can be integrated 
inwards from a point very near the surface, using the above 
expansion. By adjusting £ s the solution that has the correct 
asymptotic behaviour as £ — + can be identified. By this 
method satisfactory solutions have been obtained without 
internal critical points. 

Examples of these solutions are shown in Figs 1 and 2. 
There are several points of interest. For 7 < 5/3, the solution 
matches only asymptotically on to that of NY, deviating 
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Figure 1. Time-dependent self-similar solution for 7 = 4/3 and a = 0.1. The surface occurs at £ s sa 0.5101. The quantities e* and 
F e „, defined by e = e,(GM) 2 / 3 i- 2 / 3 and F s = F £t ,Jt~ 2 , are dimensionless versions of the specific total energy and the total energy 
flux. Note that £ 3 / 2 £7» is the ratio of the angular velocity to the Keplerian angular velocity, while £e* is proportional to the ratio of the 
temperature to the virial temperature. The dashed lines represent the steady solution of NY, which is the limiting form of the present 
solution as £ — > 0. 
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Figure 3. The variation with time of the ratio of the angular 
velocity to the Keplerian angular velocity at any fixed radius, for 
the same solution as in Fig. 2. 

substantially from it at larger radii. In the case 7 = 5/3 
the flow is differentially rotating, although it is highly sub- 
Keplerian at small radii. For a <C 1, as in these examples, the 
radial velocity is everywhere highly subsonic with respect to 
a self-similarly expanding coordinate system. Although the 
temperature is comparable to the virial temperature over 
much of the flow, this does not drive a supersonic outflow at 
large radii. Finally, the flows are bound to the central object 
in the sense that the total energy is negative. As in a thin 
disc, there is a net outflow of energy from the centre of the 
disc because the central object is accreting bound matter. 



4 DISCUSSION 

In this paper the equations of time-dependent quasi- 
spherical accretion have been solved by semi-analytical sim- 
ilarity methods. The solution describes asymptotically the 
fate of an initially cool and narrow ring of viscous fluid or- 
biting around a central mass when it is unable to radiate 
efficiently. The flow is finite in size, mass, energy and an- 
gular momentum and avoids many of the strictures of the 
steady self-similar solution of NY. 

In particular, a solution was found for the important 
case 7 = 5/3 that has differential rotation and viscous dissi- 
pation, and resembles a disc at large radii. No barrier to ac- 
cretion was found to exist in this case. At a fixed radius, the 
flow gradually becomes increasingly sub-Keplerian (Fig. 3). 
The limit 7 — > 5/3 was shown to be singular so that the 
Bondi solution is not attained. This conclusion applies also 
to steady solutions that are not self-similar. 

The flows are bound to the central object in the sense 
that the total energy is negative. Indeed, for 7 > 4/3 the spe- 
cific total energy is negative at all radii. It is suggested that 
previous interpretations based on the sign of the Bernoulli 
parameter may be insecure because Bernoulli's theorem does 



not hold for a viscous fluid (nor for a time-dependent one, 
as here). The meaning of the Bernoulli parameter is espe- 
cially unclear in a flow that is driven specifically by viscous 
stresses. 

The present model admits, at least in principle, the pos- 
sibility of a supersonic outflow at large radii, but this was 
not found to be necessary. This suggests, contrary to Bland- 
ford & Begelman (1999), that quasi-spherical accretion can 
take place without the need for powerful outflows. However, 
because the model is topologically restricted, this conclusion 
must be regarded as tentative and one cannot exclude the 
possibility that an outflow would nevertheless occur from the 
hot inner part of the flow if the one-dimensional assumption 
were relaxed. This might be answered by restoring the lat- 
itudinal dependence in the equations and finding whether 
the analogous time-dependent self-similar solution includes 
a bipolar outflow. However, this would require the solution 
of a non-linear free-boundary problem of elliptic (or possibly 
mixed) type, which is technically demanding. It is perhaps 
more likely that the initial-value problem with latitudinal 
dependence could be solved using a time-dependent numer- 
ical method. Indeed, Igumenshchev, Chen & Abramowicz 
(1996) have already performed a calculation of this type, 
adopting 7 = 4/3 and starting from a thick torus. They 
found rather complex behaviour, including convection, but 
no evidence for bipolar outflows. 

The time-dependence of the flow is not essential to the 
conclusions of this paper. However, the present method of 
investigation allows a description of generic properties of 
solutions of the initial-value problem that is essentially free 
from infinities and contains no undetermined parameters, 
only the physical ones 7 and a. Questions concerning the 
boundness of the flow and the necessity of outflows can be 
addressed with greater confidence within such a framework. 
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